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ABSTRACT 

In the framework of the Tellus Project of the European Com- 
munities, an algorithm has been worked out permitting the eva- 
luation of the thermal inertia and of the surface relative humidity 
of bare and scarcely vegetated soils from remotely sensed day 
and night surface temperatures. The cumulative daily evapora- 
tion can also be calculated by this method, 

The algorithm was tested on existing data from a field expe- 
riment and was applied in a flight experiment. While the calcu- 
lated soil moisture agreed satisfactorily with the average mois- 
ture content of the top soil layer, the estimated daily evaporation 
deviated from the measured values due to simplifications inherent 
in the model. 


1, INTRODUCTION 

The countries of the European Communities carry ou. a joint investigation on the applica- 
bility of HCMM-satellite thermal imagery for soil moisture and heat budget evaluation in 
zones of agricultural and environmental interest; the HCMM-Tellus Project. In this frame- 
work, a new algorithm has been developed to convert remotely sensed day and night tempera- 
tures of bare and scarcely vegetatod areas into parameters relevant to soil moiBture. 

Such an interpretation routine calls for a process, which is the reverse of temperature 
simulation. For this purpose, and mainly in the U.K. , methods have been worked out to cal- 
culate the thermal inertia of the soil from the difference between the day and the night tempe- 
rature to be obtained by means of HCMM fl, 2, 3, 4, 5 ], We have discussed this approach £6^ 
which 1 b limited to cases of known or negligibly smal' evaporation. 

It must be noted, however, that in the above approach the available information is not 
used completely. Only one data point i.e. the day-night temperature difference is used, 
while two data, the day and the night temperature, are available. Therefore, in principle, 
two unknowns could be determined. This has been achieved in the "TEDD-US" algorithm des- 
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cribed in the present paper. This algorithm determines the thermal inertia and the surface 
relative humidity. Moreover, the total daily evaporation is obtained. 


2. TUN ALGORITHM 

The algorithm consists of a main program and of a sub-program simulating the daily 
course ol the surface temperature (Fig. 1). The main program compares cite simulated day 
and night tomperatures with the temperatures obtained from Lhc satellite. In one version it 
adjusts the values of the thermal inertia and the surf, o relative humidity until agreement 
between the simulated and the measured temperatures has been obtained. 

A second version of the main program produces a look-up table, consisting of tbu day tempe- 
ratures, the night temperatures and the dally evaporation totals, for various combinations of 
thermal inertia and surface relativ humidity. 


The simulation sub-program selves numerically the finite difference equation which des- 
cribes transient boat conduction in a homogeneous soil: 


d T _ _A d 2 T 

dt ° C .2 
da 


(1) 


(For the meaning of the symbols wo refer to the table at the end of this paper) 

As a lower boundary, a level in the soil at a depth outeide the /uach of the daily temperature 
cycle is taken. For this level the temperature has to be specified. As an upper boundary a 
level ^ii.the atmosphere is taken, whore windspeed, air temperature and air specific humidity 
u>*e measured. The atmosphere and the soil regime are coupled by a continuity requirement; 
tiic surface heat balance: 


I + G + H + LB 0 
n 

(*> 

I n = (l*a)I + I f + £1, 

w 

G = A (~) 

'dz 'z=0 

(*i) 

T. -T 

I I = 0 c — — - 

V 

(5) 

S -h.S 


LB sOL — 

r 

(6) 


The atmospheric resistance r depends on the surface aerodynamic roughness, on the wind- 
speed and on the stability of the atmosphere. It is calculated from the integrated Businger- 
Dyer profile relationships C7 J . 

The simulation sub-program may contain only two unknowns: the surface relative humidi- 
ty and the thermal Inertia. This implies that we have to accept certain physical simplifica- 
tions. The soil is assumed to be homogeneous and with a constant water content. The surfacu 
relative humidity is also assumed to be constant. Accordingly, the values of thermal inertia 
and surface relative humidity that are produced by the algorithm, should he taken as "effec- 
tive" values. 

One might question why the thermal inertia did not show up in equations (1) - (6). The 
thermal inertia, defined as: 

P * (AC) 1 
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appears as the soil property that affects the magnitude of changes in surface temperature if 
differential equation (1) is ana lyti tally solved, for simplified boundary conditions and heat 
balance. Therefore, the observod day and night temperatures depend on the thermal inertia, 
rather than on the thermal conductivity. 


In the present, somewhat more realistic description, oar system of equations (1) - (6) 
included, had to be solved numerically by liio explicit finite difference scheme of Dufort and 
Frankul, applied to a grid which expands with depth. This scheme is used to predict succes- 
sive temperature profiles. The surface temperature is iteratively solved from the surface 
heat balance equation after each prediction, 
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Table 1 reviews the essential input data, and their determination. The data have boon 
split into local data and regional data. The regional data are to bo obtained from the nearest 
weather station. The local data have to be specified for every ground plot to be interpreted. 
Of those, the surface temperatures and the albedo may be obtained by means of the HCMM 
satellite or an airborne scanner. The most cumbersome parameter to be determined s the 
surface aerodynamic roughness. It can be estimated with existing techniques from the dimen- 
sions of the actual roughness elements employing land-use maps, which may bo derived from 
airborne, or oven LANDSAT imagery, If the area is not horizontal the slope direction and 
slope dip have to be estimated from topographical maps. 


, , '1. RELATION TO "SOIL MOISTURE" 

In terms of soil moisture, the meaning of the thermal inertia and of the surface relative 
humidity of the soil is not self-evident. While the surface relative humidity duponds uniquely 
on the soil moisture potential at the surface, thermal inertia is a function of the average water 
content in the top soil. This relation, however, is not unique as it depends also on the soil 
mineral composition. 


S. TEST OF THE ALGORITHM 

The TELL-US algorithm has been tested on data from a field experiment of the USWCL at 
Phoenix, Arizona. On March 3, 1971, a test plot of hare Avondale loam was irrigated with 
10 cm of water, and loft to dry. Wo have used the data measured 2, 6 and 14 days after irri- 
gation representing subsequent stages of drying. The soil temperature measured at "1mm 
depth" was taken to be a good approximation of the surface temperature. 

The look-up table produced by the computer program for 2 days after irrigation is plotted 
as a look-up graph in Figs. 2a and 2b, In Fig. 2a the simulated day-night temperature pairs 
are plotted for the various combinations of thermal inertia and surface relative humidity, and 
connected by lines of equal thermal inertia (dashed) and equal surface relative humidity. In 
Fig. 2b, the total daily evaporation is plotted versus the surface relative humidity. Similar 
look-up graphs were produced for the other two days of the experiment. 

The measured day-night temperature pairs are plotted as crosses in Fig. 2a and the cor- 
responding values of the thermal inertia and surface relative humidity are found by interpola- 
tion. Subsequently, the total dally evaporation values are found from the graph 2b. Instead of 
carrying out this graphical interpretation, the same result may be obtained directly with the 
automatic search version of the algorithm. 

The thermal inertia found on the look-up graphs was translated into soil moisture using 
the thermal inertia-water content relation shown in Fig. 4. This relation was calculated from 
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the mineral composition of the soil. 

The noil moisture profiles measured gravimotrlcally 2 and 1*1 days after irrigation, , in- 
plotted in Figs, 3a and b, both for the day and night measurement. The water content tabula- 
ted by the TELL-US algorithm is indicated by an arrow. Thu agreement between measured 
and calculated soil moisture is excellent 2 days after irrigation (Fig, 3a), 1*1 days alter irrtgu 
tion (Fig, 3b) the calculated noil moisture represents only an "average 11 moisture content of 
the top 8 cm of the soil. The data of the 6th day, not bIiowii in this paper, have an interme- 
diate position. 

The values of the total daily evaporation as measured by a weighing ly sunnier and deter- 
mined with the TELI.-US algorithm are compared In Table II, This comparison is not yet 
satisfactory as, except for thu initial drying stage, the algorithm lends to underestimate eva- 
poration, The reason for tills tendency lies in the assumption of a constant surface relative 
humidity which causes excessive condensation to be simulated for the night. 

It may be hoped that aftor a further Improvement of the algorithm calculated evaporation 
will come closer to valuoa measured with existing field methods. 


6 . JOINT FLIGHT EXPERIMENT 

As part of the pre-launch activities of the Tellus Project, a joint flight experiment was 
carried out in mid -September 1977 on the Grondon Underwood Experimental Catchment north- 
west of London, U. K. The application of the TELL-US algorithm was one of the amis of this 
experiment. At the lime of the night-minimum and of the day -maximum surface temperature’, 
a strip 8 km long and 1.6 km wide was scanned from 1000 m altitude with m airborne 
DS-1250 multispoctral scanner having a thermal hand botween 8 and 14/un. The night was 
clear, with no wind and the day was sunny, with a light wind and occasional high altitude 
clouds, 

Simultaneously with the overflight, ground surface temperatures wore measured with e 
PRT-5 radiation thermometer (9, 5-11, 5 , 01111 , 2° FOV) on a 6 ha ploughed field of Denehwcrih 
clay. Sixty moisture samples of the top 3 cm of soil were taken in a pattern 01 four lines 
across the field, with an additional 10 core samples down to 15 cm depth, Input data lor the 
TELL-US algorithm were obtained from instrumentation installed on the field and from an 
automatic weather station at a distance of 1, 5 km. 

First preliminary re Us of the experiment are presented in Tables III and IV. Tlu* field 
war divided into two parts of about equal size, according to the two distinct temperature clas- 
ses shown by the daytime scanner image. 

The day and night surface temperatures are given in Table 111, The airborne surface 
temperatures were corrected for the effect of the atmosphere using temperature, humidity 
and pressure data from a tethered balloon 40 km from the site. The night-time correction is 
less than the day-time correction and has the opposite sign due to a temperature inversion in 
the first 100 m of the atmosphere at night. While the night temperatures measured by the 
scanner and or the ground lie in approximately the same range, day-time temperatures mea- 
sured from the aircraft exceed the ground surface temperatures by a few degrees. This dis- 
crepancy may bo duo to a difference in the temperature distribution on the rough soil surface 
as seen by the small-angle radiation thermometer on the ground and by the airborne scanner 
with a resolution element of 2. 5 x 2, 5 m. 

As shown in Table IV, the soil moisture percentage in the 0-3 cm layer seems to be 
slightly higher in the SE-part of the field, while the 0-15 cm moisture percentage is higher in 
the NW-part. The thermal inertia and the relative surface humidity in Table IV have been 
found by interpolation in the corresponding look-up table using the. ground measured surface 
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temperatures. At thin early stage o£ data evaluation it can only be observed that both the 
thermal inertia and the surface relative humidity evolve in the right direction, I. c, thermal 
inertia increases with the mean moisture lontont of the built top soil (o-lb cm), whereas sur- 
face relative humidity increases with the moisture content in the very topmost layer (0-3 cm). 


7. CONCLUSION AND OUTLOOK 

It could be shown that by moans of the TKLL-IJS algorithm tile top-soil water content, 
the surface relative humidity and the total daily evaporation may he determined from one day- 
and one night-temperature and standard ground based data. Further experimental work is 
still needed to test the complete potential of the algorithm. Part of this work is already in 
progress in the framework of the Tcllus Project of the European Communities. 

A similar or generalized algorithm might he constructed to cover vegetated surfaces. 
Those algorithms create the attractive prospect of a combined surface temperature, soil 
moisture and evaporation mapping system on the basis of satellite or airborne thermal ima- 
gery. Such data arc of high relevance in agriculture, ecology and hydrology. They are useful 
for irrigation scheduling, prediction of yields, and largo scale management of rangelands. 
They will allow agronomists and ecologists, particularly in arid and semi-arid zones, to 
follow, analyze and predict processes such as seed germination, early crop growth, drought 
damage, erosion and desert formation, Such a mapping systom could also be used for early 
warning and fighting of pests, such as the desert locust |‘8,9’J. 
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DEFINFl ION OF THE SYMBOLS 


C volumetric heat capacity of oil [j/m 3 I<] 

G soil heat flux at the surface [J/m^s] 

U atmospheric heat flux [j/m^B’J 

Ij long wave irradiance [j/m^s] 

long wavo omittance [j/m^s] 

I sliort wavo irradiance [j/m^s] 

I not radiation [j/m 2 s] 

I, heat of evaporation [j/kg] 

JL.E latent atmospheric heat flux [j/m^s] 

T (soil) temperature [l<] 

T air temperature [k] 

a 

T surface temperature [K] 

a alhedo [ fraction] 

c heat capacity of air [j/kgK.] 

h surface relative humidity [fraction] 

p soi 1 thermal inertia [j/m^s K] 

r atmospheric resistance [s/m] 

6^ air specific vapour density [dimensionless] 

S' specific vapour density at saturation, at the soil surface 

[ dimensionless] 

t time [s] 

z depth Cm] 

A soil thermal conductivity [j/rnKs] 

Q air density [kg/m^J 

£ long wave emissivity Cfraction] 



C> 





in pkoduuhjiuiy (>!<■ Tills; 
OH/GINAr, PACE IS Poult 



FIG, 1 SUBROUTINE RELATIONAL DIAGRAM 

INFORM i general information 

INPUT : input 

IN IT ! initialization 

EXEC : exocution anti manipulation of the simulation subprogram, 
automatic searcdi for the thermal inertia and surface 
relative humidity 
OUTPUT ; output 

SIMULA ! simulation sub-program 

SINIT : initialization 

DATA : weather data interpolation 

PROF : temperature profile prediction 

RAD : calculation of irradiation terms 

ARES : calculation of the atmosphere resistance 

BUDGET : surface heat budget evaluation 

STEM : iterative solution of the surface temperature 

PRSD ; intermediate output facility 

PREP : determination of day temperature, night temperature 
and total daily evaporation. 
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TABLE I. INPUT TO "TELL-US" ALGORITHM 



Datum 

How to determine 

Remarks 

1 . 

Day surface temperature 

HCMM/alrborne scanning 

Time ol observation variabl 

2. 

Night " " 

il it 

il 

11 II It It 

3, 

Albedo 

II II 

II 


d. 

Surface aerodynamic 

Estimation from land-use 

If not available, to be 


roughness 

map 


. derived from LANDSA'I 

5. 

Slo|)u diruction 

Topographical map 

imagery 

(>. 

Slope dip 

II 

II 


7. 

Sub-surfaee temperature 

Measuremont/Estimation 

At 1 m depth 

H. 

Initial surface " 

if 

II 

At midnight 

9, 

Solar irradianco 

Radiometer 

* 


10. 

Wind speed 

Anemometer 

Weather 

Every hour 



> 

11. 

Air temperature 1 

Dry and wot 

Station 


12. 

Air specific humidity J 

bulb thermo- 
meter 
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1’AIU.K II. TOTAL DAILY [EVAPORATION' , im/day] 


T- ■ r— '<»*. • - ■ 


It* v 

Measured by 


lyflimoler 

> 

*» 

-5. 8-1 


-1.13 

i-l 

-- = -t -.ST* 

*0,86 


( + : condensation) 


Cali nliitml with 
TK1.L-US 

-6.-18 
-0. 59 
+ 1,0} 


IA1M.K m, JOINT FLIGHT EXPERIMENT - MEASURED HtiKFAl K 
TV M PE RAT U RES OF 15 ARK SOIL 


Kiold 


Radiative Surface Ti 

mperutuiv *\ 

% • 

... 

Surface 

Ai rho rue 

» 

Airborne currei ted 

NW 

Night 

Day 

4.0 +0.6* 
24.2 +1.5 

3.0 - 4. 5 
26 - 29 

2.9 * *». 1 
27 - 50 

SK 

Night 

Day 

2.-1 +0.7 
£3,9 + 3.2 

3.0 - *1, 5 
2 5 - 26 

2,9 * 4.-1 
2-1 . 27 


( i : j; standard deviation) 



i'ALI.K IV. JOINT FLIGHT EXPERIMENT - MEASURED SOIL MOISTURE PKIU.I-.N- 
IAGL AND CALCULATED PARAMETERS RELATED 'J’O SOIL MOIH'I URK 


Field 


Soil Moisture 

Percentage (weight) 

Thermal 

Relative 



0*3 cm 

0-15 cm 

Inertia 

humidity 

NW 


9. 0+0. 9* 

32.2 + 0.7 

1700 

0. 27 

SK 


10. 1)4 1,2 

28.6 + i. 8 

1300 

0. 32 


(i : + standard deviation) 
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